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Abstract 

We present a phenomenological study of the xp dependence of quarkonium 
production in high energy proton- nucleus colhsions. The xp dependence of 
comover contributions is introduced to account for the observed quarkonium 
suppression at low xp. Combining comover contributions, nuclear shadowing 
effect, energy loss mechanism and nuclear absorption together we reproduce 
the overall xp dependence of E772/E789 data. 
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I. INTRODUCTION 



A strong suppression of charmonium (J/-?/', 'ip') and bottomonium (T15, T25+3S') produc- 
tion has been observed in p — A collisions on heavy relative to light nuclear targets flj-Q . 
Special attentions have been paid to the significant suppressions at low (small or negative) 
and high xp regions. The suppression at low Xi? is usually believed to come from nuclear 
absorption and comover contributions. The maximum contribution of nuclear absorption to 
the J/ if) suppression was studied in ref. 0, and it was found that nuclear absorption alone 
cannot account for the J/ ip suppression measured in either p — A or A — B collisions. Since 
at E772 energy most of the physical bound states are formed well outside target nucleus due 
to the Lorentz dilation of the formation time, nuclear absorption can contribute little to the 
observed suppression at low xp region. While comovers can continue to interact well outside 
of the nuclear volume, they may in principle play a more important role in proton-nucleus 
collisions, although the number of comovers in p — A reactions is not so large as in A — i? 
collisions. In the literature comover contributions are generally estimated in the central 
rapidity region without xp dependence. To our knowledge there is no direct explanation 
of the E772/E789 data at ~ and below PJ^. In this work we use a simple model to 
extrapolate the rapidity distribution of comovers for p ~ A reactions from p — p data , 
introducing the xp dependence of comover contributions, to explain the E772/E789 data. 

On the other hand, several mechanisms have been proposed to account for the suppres- 
sion at large xp . Among these, nuclear shadowing effect alone is argued to be unable to 
explain the data [^f. The intrinsic charm model was expected to explain this xp de- 
pendence. However, E789 data of J/ip production at very large xp show no evidence 



of intrinsic charm contribution. The energy loss mechanism could account for the ob- 
served suppression at large xp . Nevertheless, the amount of energy loss needed to explain 
the E772 J/ ip data appears to be significantly larger than the value determined from other 
considerations [|ll|]. In our work, we combine energy loss mechanism together with nuclear 
shadowing effect, nuclear absorption and comover contributions, and reproduce the large 
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xf data with a smaller amount of energy loss. 



II. COMOVER CONTRIBUTIONS 

In the case of quarkonium production, e.g., J /ip production, in proton- nucleus collisions, 
comoving secondary particles can scatter with the cc pair, and contribute to the observed 
J/ip suppression. These comover contributions have been studied in the literature |T^-[l4|,^. 
However, since the shape of the rapidity density of the comovers in inclusive J/ ip production 
for p-A collisions is unknown, comover density is usually estimated in the central rapidity 
region, and taken as a constant independent of the comover rapidity. In this work we will 
use a simple model to extrapolate the rapidity distribution of comovers for p-A reactions 
from p-p data 0. Assuming that comovers with rapidity close to the rapidity of a cc pair 
can cause the breakup of the cc pair, we translate the rapidity distribution of comovers 
into xp distribution, and include the resulting xp dependence of comover contributions to 



explain the E772/E789 low xp data. We will follow the work of refs. P,p!4|, and improve the 
treatment of comover contributions. 

Including only the comover contributions, the Xp dependent cross section of J/ip pro- 
duction 'm.p — A collisions can be expressed as 



dxF dxi 



J d% J dzpA^b, z)exp 1^— J dT(TcoVn{T,b)^ , (1) 



where pA is the nuclear density profile, cTco is the (cc)-comover absorption cross section, tq 
is the formation time of comovers, Tf is the effective proper time over which the comovers 
can interact with cc pair, v is the relative velocity of cc with the comovers, and n(r, b) is the 
density of comovers at the proper time r and impact parameter b. 

Assuming Bjorken's hydrodynamics |15|, the comover density varies with the proper time 



according to. 



n(r) = (2) 



3 



The comovers consist mostly of vr, p and uo mesons. We relate n, the density of comovers, 
to rich, the density of produced charged particles, hy n = frich- For example, / ~ 1.5 if 
comovers are all tt's with equal numbers of 7r+, tt^ and vr". 

The charged particle density at tq can be related to rapidity distribution by 

nch{To,b) = f^, 3 

where (Tj„ ~ 30 mb is the nucleon-nucleon inelastic scattering cross section. 

Now we use the so called independent cluster model [|l6l to extrapolate the charged 



particle rapidity distribution for p — A reactions from p — p data . 

We fit the data at ^/s = 45.2 GeV, which is close to the E772 energy, in the form. 



Wexp [-{r] + riof/5] + Wexp {-(r/ - riof/S] (4) 



drj 

with W=1.75, r/o=l-5 and 5=2.9. 

In the framework of the independent cluster model, we have 

^^^^^^ = i^{b)Wexp [-{7] + + Wexp {-(r/ - , (5) 

where uib) = ain J dzpA{b,z) is the number of target participants at the impact parameter 
b. 

In order to introduce the xp dependence of comover contributions, we assume that those 
comovers with rapidity close to the rapidity oi a. J/ip can cause the breakup of the J/V', 
and we relate the rapidity of J/ip in the center-of-mass system to xp-, which is defined as 

IF = Pil(n),na., by 

^ 2 \E* -Plj' ^ ' 

where E* and P£ are the energy and longitudinal momentum of J/ip in the c.m.s.. 

In our approach to the xp dependent comover density two main parameters are / and 
tq. In Fig.l we show the xp dependence of comover densities for various sets of / and tq. 
Note that the dotted curve with / = 1.5 x 0.6 and tq = 0.8/m is plotted in view of the 



argument of ref. [0, which represents a comover density too low to explain the data. In our 
further calculations we will fix / = 1.5. The comover formation time has been estimated to 
be about 0.4 — 1.2 /m/c, the time scale of both quark and hadron formation in the early 
stages of the collision p!7|JT8| . It can be seen from Fig.l that the comover density is sensitive 
to the choice of tq. In the next section we will vary tq and cXco to attribute the observed 
suppression at low xp region to comover contributions. 

Since comovers can interact with the produced heavy quarks Q or Q before the quarko- 
nium bound state can be formed, and the interaction is independent of whether or not the 
heavy quark pair is produced as a color octet or singlet, we expect that for charmonium 
production the (cc)-comover interaction cross section is the same for the J/ip and ip' bound 
states. Furthermore, the comover contributions are expected to be independent of the fi- 
nal hadronic size. As a result, although we expect the occurrence of significant comover 
absorptions with strong nulear dependence, the J/ip and ip' production would have the 
same nuclear dependence under comover absorption, which is consistent with the earlier 
experimental observation H|. 



III. COMBINATION OF MECHANISMS, NUMERICAL RESULTS AND 

DISCUSSIONS 

In this section we combine comover contributions, nuclear absorption, nuclear shadowing 
effect and energy loss mechanism to reproduce the overall xp dependence of E772/E789 
data on J/ ip and T production. 

We will calculate xp dependence of a for J/ip and T production in the parametrization, 

dxp dxp 

For J/ip production, we make use of the fact that about 60% of the final J/ip^s are produced 
directly, about 30% come from Xc intermediate states and the remaining 10% are produced 
through the decay of ip'- Analogously, about 41% of the final Ti5's are produced directly, 



5 



about 35% and 18% come from the decay of Xb and Xb^ about 5% and 1% are produced 
through the decay of T25 and T35 |]I9| . 

The cross section of quarkonium production in p-A colhsions is 



J d% J dzpA{b,z)exp i- / dz' pA{b, z')aabs{z' - z) 

xexpj— y dTacoVnij.h)^ (8) 

where several mechanisms have been included. 

The averaged nucleon-nucleon cross section oip — A reactions, aPP* in eq. (^), is the same 
as the bare nucleon-nucleon cross section except that in aPP' the distribution functions of the 
gluon, quark and antiquark within the target nucleus are nuclear modified. We follow the 
work of ref . to incorporate the shadowing effect into aPP' , and use parameters determined 
in ref. §. 

The first exponential in eq. (§) represents the contribution of nuclear absorption. We 
take the values of parameters in the following way: the formation time of charmonia and 
bottomonia, = 0.89 /m, tx = 0.76 /m, etc. , is taken from ref. the absorption cross 
section between resonances and nucleon is chosen as 

c^RN = CF^N{rR/r^f, (9) 



with fj^AT = 6 m6, and the radii of resonances taken from ref. |2T|]. In particular, we want 
to point out that with the reason discussed in the first section, the nuclear absorption can 
contribute very little suppression even with a rather large absorption cross section — 
6 m6, as indicated by the fact that the nuclear dependence of J/ip and ip' production is the 
same within errors 0. 

We follow the work of ref. [|T^ to include the contribution of energy loss mechanism 
in our calculations. In ref. |]T0[ only energy loss mechanism is taken into account, and an 
upper bound of energy loss is obtained, dE/dZ ^ 1.5 GeV/fm. Moreover, the energy loss 
contribution in ref. |]TD[ has not been normalized, which implies an integrated suppression. 



In our calculations, we complete the normalization in a similar way as in ref. ^2 



da^^ 1 da^"^ 

-J [Xp) \normal.= 1 {Xf/o), (lO) 

dxp a dxp 

where a = ^^^^^^ , and Axp is the shift of due to the energy loss. Since nuclear 
shadowing effect can also contribute to the quarkonium production suppression at large xp, 
in our calculation a value of dE/dZ ^ 0.5 GeV/fm is found to be able to fit the data quite 
well, which is consistent with the value obtained from other considerations [|TT[] . 

The comover contributions are included in the second exponential in eq.( |^). We take 
the effective proper time tj = vq/cs, where tq is the projectile radius which we take to be 
1.2 fm, and c, ~ 1/^3 

The results of our calculations are compared to E772/E789 data for J/ip production at 
800 GeV/c in Figs. 2 and 3, with different sets of dco and tq. We find that (Xco = 4.0 mb and 
tq = 0.8 fm/c is the best fit, which corresponds to a comover density with maximal value 
about 6 fm~^. Fig.2 shows that a lower density (larger tq) is not sufficient to explain the 
data. In Fig. 3 we see that a comover absorption cross section aco smaller than 4.0 mb is also 
insufficient to fit the data. Note that the dot-dashed curve in Fig. 3 corresponds to the case 
without comover contributions. The curve illustrates a very small contribution of nulear 
absorption at low xp, as we expected. One can find in Fig. 3 that although the number 
of comovers is not very large in p — A reactions, comover interactions have a much more 
important effect than nuclear absorption. 

A comparison of our calculations for T production with the E772 data is shown in 
Fig. 4. The dashed curve with aco = 4.0 mb and tq = 0.8 fm/c is somehow lower than the 
experimental points at xp — 0.3. This may imply our over-estimate of comover contributions 
when we fit the data at low xp region. 

In the literature studies of nucleus- nucleus reactions indicated that the comover density 
at an early stage of the collision is very high, possibly from 1 to 5 fm~^. In our work, a 
comover density with a maximum~ 6 fm~^ is used when we fit the E772/E789 data at low 
Xp. With such a high density one would rather think of these comovers as light quarks and 



antiquarks than as pions and low-mass resonances [Q. However, whether it is possible to 



obtain such a high density in proton-nucleus coUisions appears questionable. Presumably 
some new mechanism might also contribute at low xp region. 

According to ref. ||23[, absorption in confined hadronic matter is excluded as a possible 
cause of the J/ip suppression observed in nucleus- nucleus collisions. Since at low xp region 
neither energy loss nor nuclear shadowing effect is important, it seems difficult to under- 
stand the low xp suppression in proton- nucleus collisions, in which no deconfined matter 
is expected. Our work shows that comover interactions play an important role at the low 
Xp region, and a comover absorption cross section as large as 4 mb is needed to confront 
the E772/E789 data. In addition, a recent work observed that nuclear antishadowing leads 
to enhanced J/ip production at low xp, which is in contradiction with the E772/E789 data 



2J]. These issues suggest that the low xp suppression in p — A reactions deserves careful 
studies, and possibly requires some new mechanism. 

To summarize, we have reproduced the overall xp dependence of quarkonium produc- 
tion in proton-nucleus collisions, by combining comover contributions, nuclear shadowing 
effect, nuclear absorption and energy loss mechanism together. In particular, by including 
the Xp dependence of comover distribution we show that comover interactions have a very 
important effect when one attempts to explain the observed quarkonium suppression, es- 
pecially at low Xp region, in proton- nucleus collisions. However, our work indicates that 
a fairly high density of comovers is needed to explain the E772/E789 low xp suppression, 
which might imply some new mechanism occurs at this region as well. 
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Figure Captions 



Fig.l. xp dependence of comover densities for different sets of / and Tq. 

Fig. 2. xp dependence for J/ifj production from E772/E789 [3,4]. The curves compare 
our calculations to the data, with a fixed aco = 4.0 mb and various values of tq. 

Fig. 3. xf dependence for J/t/j production from E772/E789 [3,4]. The curves compare 
our calculations to the data, with a fixed tq = 0.8 fm/c and aco = 4.0 mb, 2.0 mb. The 
dot-dashed line corresponds to the case without comover contributions. 

Fig.4. Our calculations oi xp dependence for bottomonium production with different 
sets of (7co and tq are compared to the E772 data [3] . 
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